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Abstract 

The ground state energy and entropy of the dilute mean field Ising model is computed exactly 
by a single order parameter. An analogous exact solution is obtained in presence of a magnetic 
field with random locations. Results allow for a complete understanding of the geography of the 
associated random graph. In particular we give the size of the giant component (continent) and 
the number of isolated clusters of connected spins of all given size (islands) . We also compute the 
average number of bonds per spin in the continent and in the islands. Then, we tackle the problem 
of computing the free energy of the dilute Ising model at strictly positive temperature. We are 
able to find out the exact solution in the paramagnetic region and exactly determine the phase 
transition line. In the ferromagnetic region we provide a solution in terms of an expansion with 
respect to a second parameter which can be made as accurate as necessary. All results are reached 
in the replica frame by a strategy which is not based on multi-overlaps. 

PACS numbers: 05.50.+q, 64.60.De, 75.10.Hk, 81.05.Kf 
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1 - INTRODUCTION 



The study of the dilute mean field Ising model can be tackled using replica approach 
which has been successfully applied to many disordered systems, the most celebrated being 



the SK model 



2[. Nevertheless, the method encounters serious difficulties when applied 



to many other systems because of a proliferation of replica order parameters as, for example, 



the multi-overlaps in dilute models [3 



In this paper we use a strategy proposed by Monasson which we already applied 
to the dilute Ising model for which we computed the exact entropy and energy in the 
vanishing temperature case. The Monasson strategy is very general and shows how <Z n > 
(Z is the partition function of a generic spin model with disorder) can be expressed in 
terms of a maximum with respect to the possible values of 2 n positive densities x(a) (where 
a = (cr 1 , a 2 , , a n )) with the constraint J2a x ( a ) — 1- 

In next section, following we derive the free energy of the dilute Ising model in terms 
of the densities x(a) for a generic temperature and also in presence of an external magnetic 
field with random locations. In section 3, we consider the zero temperature case, in absence 
of magnetic field and we show that the model is formally equivalent to a 2 n Potts model. This 
observation allows to express the x(a) in terms of a single order parameter and the exact 
solution is easily obtained. The results C01 „ C1 de with those found in fl where the approach 
is not based on replicas. We also show that the exactly computed internal energy and 
entropy allows to draw conclusions about the geography of the associated random graph. In 
particular we determine the size of the giant component (continent), the average size of the 
isolated clusters of connected spins (islands) and the average number of bonds per spin in the 
continent and in the islands. In section 4 we compute exactly the ground state energy and 
entropy in presence of a randomly located magnetic field. This allows for a more complete 
description of the geography, in particular, the number of isolated clusters of connected spins 
of all given size (islands) is determined. In section 5 we extend the scope by considering the 
same model in the case of a strictly positive temperature, where a second order parameter 
is needed. We are able to find out the exact solution in the paramagnetic region and exactly 
determine the phase transition line. In the ferromagnetic region we provide a solution in 
terms of an expansion with respect to the second order parameter which can be made as 
accurate as necessary. Conclusions and outlook are in the final section. 
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2 - THE DILUTE FERROMAGNET AND THE REPLICA SOLUTION 



The dilute ferromagnetic Ising model is characterized by a number M of order N of non 
vanishing (and positive) bonds connecting pairs of spins. The partition function can be 
written as 

Z = X^exp 1/3 Yl K ij°i°j (!) 

# V »>j / 

where the sum J2# goes on the 2 N realizations of the N spin variables and the are 
quenched variables which take the value 1 with probability and otherwise. The dilution 
coefficient 7 may take any positive value and the number M of bonds is ± o(y/~N). This 
is at variance with the fully connected model where the number of non vanishing bonds is 
N(N — l)/2 and it is the reason why the dilute ferromagnetic has a richer behavior. 

An alternative definition of the partition function can be obtained by taking a number 
of bonds exactly equal to M where M can be both a deterministic number proportional 
to iV or a quenched random variable with average proportional to N and fluctuations of 
order yN. In this context the pair of spin to be connected to a given bond is chosen 
randomly (4,8]. All these models and model (pQ) are thermodynamically equivalent, i.e. the 



intensive thermodynamical quantities are the same. 

Models with lesser dilution 9, [uj] (number of bonds of order N e with l<e<2 ) have 
been also considered, but in this case the thermodynamics is the same of the fully connected 
model. 

The problem of evaluating (CQ) via the replica method, can be approached in many different 
ways. A possibility is to follow the same path which has permitted the solution of the SK 
model using multi-overlaps as r epl,ca order parameters. This strategy leads to a proliferation 
of parameters and it is unable to give exact answers |5|. In [7], we considered a different 
strategy, which is based, instead, on the densities x(a) which are 2 n parameters since a = 
(o" 1 ,^ 2 , , <7 n ) may assume 2 n possible values. 

In [7j we computed the partition function and, therefore, the free energy F obtaining 

-/3F=lim^ + ^log(cosh(/3)) (2) 

where \l/ n is given by 
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tt n = max [ -1 + | ]T x(a)x(r) + tanh(/3)a a r a ) - £ log(x(«7)) ] (3) 

where the sum Y^a r g° es on the 2 n possible values of the variable cr and the the 2 n possible 
values of the variable r. The maximum is taken with respect the 2 n densities x(cr) with 
the constraints < x(cr) < 1 and J2a x { ar ) = 1- The highly non trivial problem is the 
maximization of the x(a) which, in principle, can be found by a proper parametrization. 

The simplest case turns out to be the high temperature paramagnetic region where the 
maximum is reached when x(a) = ^ for all a, so that 

-/3F = log(2) + |log(cosh(/3)) (4) 

We will show in section 5 that the paramagnetic region is the high temperature region 
given by 7tanh(/3) < 1 (see also [4j, 

In 7| we also considered the dilute ferromagnetic system in a magnetic field with random 
locations. The partition function of this model is 

Z = ^2 exp /3 J] Kijaiaj + (3 ^ (5) 

# \ i>3 i J 

where Kij are the previously defined quenched variables and the hi take the positive value 
h with probability 5 and otherwise. 
In this case we get 

-pF= lim + ^ log(cosh(/3)) + 5 log(cosh(/3/i)) (6) 

n->o n 2 

where \l/„(5) is obtained by maximizing an expression which is the same of that in (EJ) plus 
the extra term 

5>(a)log (l -5 + 5Yl(l + tanh{(3h)a a yj (7) 

cr \ a / 

where the sum J2<r g° es on the the 2 n possible values of the variable cr. Obviously, \I / „(0) = 

Before trying to compute *ff n (8) in the general case, we focus on the zero temperature case, 
which is exactly solvable both in the vanishing (next section) and non vanishing (section 4) 
magnetic field case. In both cases many conclusions may be drown concerning the structure 
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of the random graph, both for what concerns the giant component (continent) and the small 
islands of connected spins. 



3 - ZERO TEMPERATURE AND ZERO MAGNETIC FIELD 

Let us start by considering the case of vanishing temperature and vanishing magnetic 
field (5 = 0). In the vanishing temperature limit one has tanh(/3) = 1 and, therefore, 
Yl a (l + tanh(/3)cr Q r a ) becomes n a (l + o a T a ). This term equals 2 n when all a a equals r a 
and vanishes otherwise. Therefore, expression (J3)) becomes 

* n = max [ -1 + \T £ x{af - £ x(a) \og(x(a)) ] (8) 



which is a standard 2 n -components Potts model [U|]. The solution is known and can be 
found assuming that 2 n — 1 quantities x(a) take the value and one takes the value 
1+( ' 2 2 ,~ 1 ^ . The state with different value can be any of the possible 2 n , we assume that is 
the one with all cr a = 1 for all a. We can write: 

j(g) = l^ + grc4±o (9) 

which satisfy the constraint J2a x ( cr ) = 1- 

In [3] we used this solution in order to find the exact expression for the ground state 
internal energy E 

E = -Z (10) 



while the entropy S was found to be 

S = log(2) max [^(6 2 - 1) + (1 - 6)(1 - log(l - 6)) ] (11) 
The maximum is reached in 9 C given by the equation 

exp(- 7 # c ) = l-6 e (12) 

This equation has a single non negative solution 9 C = if 7 < 1 and one more non trivial 
positive solution if 7 > 1 which corresponds to the maximum. Therefore, at temperature, 
for 7 < 1 the system is in a paramagnetic phase (6 C = 0) while for 7 > 1 is in a ferromagnetic 
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FIG. 1: Entropy (dashed line) and order parameter 8 C (full line) as a function of the dilution 
coefficient 7 at temperature and magnetic field. The transition is at 7 = 1 where the first 
derivative of 9 C and the third derivative of the entropy are discontinuous. 



phase (6 C > 0). The parameter 8 C is the magnetization of the system, but we will see that 
it has a simple interpretation in terms of underling random graph and that its discontinuity 
corresponds to the percolation transition generated by the ferromagnetic links. 
Using equation ()12p we can rewrite the entropy (TTTT) as 



^ = log(2)[(l-I) + ( 7 -l)^-^ 6 2 ] (13) 

The entropy S and the order parameter 9 C are plotted in Fig. 1 as a function of the dilution 
coefficient 7. At the transition value 7 = 1, the first derivative of 9 C and the third derivative 
of the entropy are discontinuous. 
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Indeed, these results have an interesting interpretation in terms of random graph ge- 
ography. In fact, at zero temperature, the magnetization is entirely due to the giant set 
(continent) of spins all connected, while the small sets of a single or a few interconnected 
spins isolated from the others (islands) cannot contribute to the magnetization. Since all 
spins in the continent must be directed in the same direction in order to minimize the ground 
state energy, their magnetization must be 1 and, therefore, their number must be 9 C N. This 
also means that the total number of spins in the islands must be (1 — 9 C )N. In particular, 
when 7 < 1 there is no continent (6 C = 0) and all the spin are in islands, on the contrary, 
all spins are in the continent only when 7 — » 00 (in this limit 8 C = 1). 

It is also possible to count the number of islands, in fact, any of the islands of connected 
spins can live only in two configurations: all spins up or all spins down. Therefore, any 
island contribute by a log(2)/7V to the (intensive) entropy. The continent also contributes 
by a single log(2)/7V but this is irrelevant in the thermodynamical limit. Therefore, the 
number of island is simply log / 2 ) ^- 

Let us call a(l)N the number of islands of size / (islands made by a number / > 1 of 
interconected spins which are isolated form all the others). Then, the total number of islands 
can be written as J2i>i oi(l) N. Furthermore, since an island of size / contains by definitions 
I spins, the total number of spins in the islands can be written as J2i>i cx(l)lN . 

We have seen that the contribution of any island to the entropy is log(2)/JV, therefore, 
we can write 

S = log(2)5>(J) (14) 
Analogously, since the total number of spins in all the islands is (1 — C )N, we can write 

1-0 C = 5>(Z)Z (15) 
1>1 

We can also easily compute the average size <l> of the islands (the average number of 
spins per island) since it is given by the ratio between the total number of spins in all the 
island J2i>i a (l) IN divided by the number of islands J2i>i a (0 N. We obtain 

< l >= (1 -T g < 2 > (16) 
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FIG. 2: Average size <l> of the islands (full line), average number of bonds per spin in the islands 
4>/(l — 9 c ) (dashed line) and average number of bonds per spin in the continent (7/2 — c/>)/# c (dotted 
line). The average number per spin in the continent is computed only when the continent exists 
{6 C > 0), that's why the dotted line begins in 7 = 1. 

The average size <l> of the islands is plotted in Fig. 2, notice that it has a maximum 
at the transition 7 = 1 where it equals 2. This means that the average dimension reaches a 
maximum when the continent begins to exist and then decreases because larger island are 
more easily absorbed by the continent. 

Another simple consideration allows to compute the number of bonds <pN in the islands. 
In fact, a moment of reflection shows that the number of islands forming loops is irrelevant 
in the thermodynamical limit. Therefore the spins of an island of size / must be connected 
by / — 1 bonds which implies 



S 



= £>(O (i-i) = i-e c 



s 



(17) 



log(2) 



Since the total number of bonds is ^N, we also have that number of bonds in the continent 
is (7/2 — <f>)N. Finally, we can compute the average number of bonds per spin in the islands 
as 0/(1 — C ), while the average number of bonds per spin in the continent is (7/2 — <p)/8 c . 
This two quantities are plotted in Fig. 2. The average number of bonds per spin in the 
islands follows the average dimension of the island and it has a maximum in 7 = 1 where it 
is the value 0.5. This is coherent with the fact that an island of size two as a single bond. 
The average number of bonds per spin in the continent is defined only when 7 > 1 because 
this the region where the continent exists. At 7 — > 1 it has its minimal value 1 coherently 
with the fact that the any of the spin must be connected at least to another spin, which is, 
in turn connected to a third one and so on. 

In next section we will consider again the case of vanishing temperature, but in presence of 
a magnetic field. This will allow for a better comprehension of the random graph geography. 
In fact, we will be able to compute all the a(l), for the moment we just remark that the 
number a(l)N of islands of size I (isolated spins) is iVexp(— 7). This can be computed by 
considering that exp(— 7) is the thermodynamical limit of (1 — ^) 7V_1 which, in turn, is the 
probability that a spin has no positive bonds connecting to all the others. 

4 - ZERO TEMPERATURE AND NON ZERO MAGNETIC FIELD 

In the vanishing temperature limit one has tanh(/3) = tanh(/3/i) = 1 and *& n (5) becomes 



*„(5) = max [-Uj2 n £ x{af+Y. x t°) ^g(l-6+6 H(l+0)-]T x(a) hg(x(a)) ] (18) 



This is again a Potts model whose solution can be found assuming that the 2 n densities 
x(a) have the form 

n 

With some simple algebra [7J we can compute the ground state internal energy E(5) 



a 



E(S) 




(19) 



while the entropy S(5) is 
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dilution coefficient 

FIG. 3: Magnetization 9 C {5) as a function of the dilution coefficient 7 at three different values of 
the magnetic field concentration: 5 = (full line), 5 = 0.1 (dashed line) and 5 = 0.2 (dotted line). 
One can easily remark the absence of transition for strictly positive values of 5. 



S(S) = log(2) max [ ?-(9 2 - 1) + (1 - 9) log(l - 8) + (1 - 0)(1 - log(l - 9)) ] (20) 

These two quantities obviously coincide with those computed in the previous section when 
5 = 0, i.e. E(0) = E and S'(O) = S. The maximum is reached in 9 C {$) given by the equation 



(l-5)exp(- 1 9 c (5)) = l-9 c (5) (21) 

At variance with the magnetic field case, this equation has a single non negative solution 
9 C {5) for any value of 7 and, therefore, the transition disappears. Also for this quantity one 
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has by definition 6 C (0) = 8 C where 9 C is the zero magnetic field magnetization computed in 
previous section. Using fl2Tj) . the entropy S(5) can be rewritten as 



S(5) = log(2)[(l - 1) + (7 - IMS) - le c (5f + (1 - 6 C (5)) log(l - 5)] (22) 

In Fig. 3 we plot the magnetization 9 C (5) as a function of the dilution coefficient 7 at 
three different values of the magnetic field concentration: 5 = 0, 5 = 0.1 and 5 = 0.2. One 
can easily remark the absence of transition for strictly positive values of 5. 

Finally we compute the entropy response a (5) to magnetic concentration 5 as 

... dS{8) dS{6) dS(S)d9 c dS{8) 

a{6) = ~dT = —d5 00-^6 = —BT (23) 

which gives 

•V = ^ (24) 
In order to use the above results to obtain more information concerning the geography 
of the random graph, we preliminarily observe that the difference between 8 C (5) and 6 C is 
only due to the effect of the magnetic field on the islands while the magnetization of the 
spins of the continents remains unchanged and equal to 1. The same can be said for what 
concerns the entropy, in fact, the contribution to the entropy coming from the continent 
remains equal to zero. 

It is easy to understand the effect of the magnetic field on a given island: if none of the 
component spins get a positive magnetic field, the island can still live in two configurations: 
all spins up or all spins down. For this reason its contribution to the entropy remains 
log(2)/iV and the contribution to the magnetization remains equal to zero. But if one or 
more spins of the island get a positive magnetic field, then all spins of the island must be 
up, which implies that the contribution to the entropy vanishes and the contribution to the 
magnetization rises from to l/N where I is the size of the island. 

The interesting point is that while magnetization 9 C N has a geometrical interpretation, 
simply representing the number of spins in the continent, the magnetization 6 C (S)N is the 
sum of the number of spins in the continent plus the number of all spins in the islands 
frozen by the magnetic field. The entropy S(5)N, on the contrary, is simply the number of 
unfrozen islands multiplied by log(iV). 
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We notice that the probability that none of the spins of an island of size I get a magnetic 
field is (1 — 5) 1 while the probability that at least one has a magnetic field is 1 — (1 — 5) 1 . 
The above arguments allow to say that the number of islands of size / which are not frozen 
is a(l)(l — 5) l N. Then, we are allowed to write 

5(tf)=log(2)E«(0(l-*) 1 (25) 
i>i 

1- 0c(5)=£«(O(l (26) 
i>i 

which coincide with the analogous quantities in previous section when S = 0. 

It is easy to verify that (125]) and (126]) are coherent with equation (124]) for the entropy 
response. Furthermore, they give much more information, in fact, expanding 1 — 9 C (8) 
for small values of 1 — 5 we have that the coefficients of the expansion are a(l) I. As a 
consequence, we can have in principle all the density numbers a (I) of islands of size I. 

In order to perform this expansion we can use iteratively equation (|2ip and we obtain 



a(l) = exp(-7) 

(27) 



a(2) = ± 7 exp(-2 7 ) 



a(3) = | 7 2 exp(-3 7 ) 
a (4) = | 7 3 exp(-4 7 ) 
while, in general, for any I > 2, one can deduce from fl2T]) the recursive formula 

«(0 = y «p(-7) ElJ^E A «&) k (28) 
1 k=i K - {k} i=i 

where the second sum goes on all the possible values of the k strictly positive integer numbers 
Zj which satisfy the constraint l\ + 1% + • • • + 1^ = I — 1. This recursive formula gives 

«(/) = ^y- 1 exp(- 7 /) (29) 

indeed, we are not able to rigorously prove ( 129]) from (128]) . nevertheless, we checked that 
it is exact for all I between 2 and 10. Much more important, we were able to numerically 
compute J2i>i a (0 an d J2i>i a (0 I an d verify that they coincide, respectively, with S and 
1 — 6 C . In particular J2i>i a (0 I equals 1 for any 7 < 1. 

Notice that the expression a(l) — exp(— 7 ) was already found in previous section using 
very simple arguments. Also notice that in the 7 — > limit all a (I) vanish except a(l) which 
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goes to 1. This is simply due to the fact that in this limit all spins are isolated which means 
that only islands of size 1 exist. 

The description of the geography of the random graph is now completed, since we are 
able, for any value of 7, to give both the dimension of the continent 6 C and the number a (I) 
of islands of any fixed size. Some relevant intuition about this geography can be obtained 
considering the islands index r{m) defined as 



which is the number of islands of size I < m divided by the total number of islands. The 
index r(m) goes to 1 for large m. In Fig. 4 we plot this quantity for m = 1, 2, 3 and 4. The 
index approaches 1 for large values of m but, for m = 4 it is still about 0.94 when 7 = 1. The 
lower values of r(m) at 7 = 1 mean that the largest islands appear around the transition. 

5 - POSITIVE TEMPERATURE 

For a positive temperature, in absence of magnetic field, the free energy is given by ([2]) 
once the proper maximum in ([3]) is found. We have seen that in the case of vanishing tem- 
perature the exact solution can be found because the problem is equivalent to a Potts model, 
on the contrary, for positive temperature, there is not an analogous equivalence. Therefore, 
we have to find the correct parametrization of the x(a) making a decision motivated by 
the physics. First of all we remark that for positive temperature, the islands cannot con- 
tribute to the magnetization, for the same reason why they do not contribute in the zero 
temperature case. In fact, being isolated clusters of spins, they can be orientated in both 
directions. So, we can already conclude that only the spins in the continent may contribute 
to the magnetization. Let us assume that the magnetization of a spin in the continent is A, 
then it is natural to assume that 



where 9 C is the continent size given by (j!2p and < A < 1 is the parameter to be max- 
imized. Notice that the total magnetization is m = 8 C X, being simply the product of the 
magnetization of a spin in the continent A and their density 9 C . In fact, 




(30) 




(31) 
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FIG. 4: Number of islands r(m) of size I < m divided by the total number of islands for m = 1, 2, 3 
and 4 (from below). The index approaches 1 for larger values of m. The lower values of r(m) at 
7 = 1 mean that the largest islands appear around the transition. 



for any of the replicas, i.e. any 1 < a < n. 

It is clear, at this point, that the really crucial ansatz is that the magnetization is a 
self- averaging quantity, otherwise we should chose an expression analogous to (l3~Tj) but with 
a different A„ correspondingly to any a a . This breaking of the symmetry would have implied 
a different magnetization m a for any possible 1 < a < n. 

With the choice ( 13 ip we have \l/ n = max^ $ n (A), where 




(32) 



(J 
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$ n (A) = ^ 6 2 [(1 + A 2 tanh(/3))" - 1] - 5>(<r) log(x(a)) (33) 

and where is given by ( 15T1) and the sum goes on all the 2 n possible realizations of a . 

In the zero temperature case (tanh(/3) = 1 ) and vanishing n, we can easily check that 
the maximum is reached in A = 1 and we recover the zero temperature solution of section 
2. This check can be performed by simply verifying that the derivative with respect to A of 
$„(A) equals when A = 1, tanh(/3) = 1 and n — > 0. 

To find the proper maximum and perform the n — > limit is much more difficult when 
the temperature is positive i.e. tanh(/3) < 1. The problem is that we were unable to find an 
analytic expression for any n of the sum Y^a x { a ) l°g(x(<7)) when < A < 1. Nevertheless, 
in order to find the behavior of the free energy in the paramagnetic region where A = 
and just around the transition where A 1 it is sufficient to expand $ n with respect to A. 
The minimum order of the expansion is the fourth, in fact, lower orders would be unable to 
provide a maximum. 

Let us rewrite x(a) = ^-(1 + y(cr)) where y(a) = c (Jla=i(^ + ^ ~ 1)> then notice that 
y(c) is a polynomial of A with no terms of order 0. Then we can expand Yla x ( a ) ^°s( x ( (7 )) 
to the fourth order in |/(cr), the approximate expression can be easily computed in a compact 
analytic form in terms of n and A and coincides with X^o- zi^) l°g( a; ( cr )) f° r all terms of order 
4 or less of the expansion in A. 

This analytic approximation can be inserted in (I3"3"|) replacing X^o-^K ") l°g(:r(cr)) and we 
obtain an expression whose expansion coincides with that of $ n (A) up to terms of order 4 
in A. At this point, if we retain only these coinciding terms and we expand to order 1 in n, 
we obtain 



$ n (A) ~ n9, 



y( 7 tanh(/3) - 1) - ^(3 7 tanh(/3) 2 - 3 + Q6 C - 26, 



+ nlog(2) (34) 



The maximum is reached in A = for 7tanh(/3) < 1 and in A = A c given by 



7 tanh(/3) - 1 



■v 2 _ 7 tmm^; - ± 

c 7 tanh(/3) 2 - 1 + 29 c - f # c 2 { } 

when 7tanh(/3) > 1. Notice that in this region 7tanh(/3) 2 > therefore the denominator 
in ( 155]) is larger then — 1 + 29 c — %Q 2 C which, in turn, is larger then 0. 
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The conclusion is that the magnetization vanishes when 7tanh(/3) < 1 and equals m c = 
6 C X C when 7tanh(/3) > 1 which implies that we have exactly determined the transition line 
7tanh(/3) = 1. 

Inserting (1551) in f[M|) we obtain 

-pf-1 + .„ g(2 , + f 7 Jff!^ + X w w 

which gives the free energy in the low (but not too low) temperature region tanh(/3) > 1. 
On the contrary, inserting A = in flM|) . we obtain 



-/3F=|log(cosh(/3)) + log(2) (37) 

which gives the exact free energy in the high temperature region tanh(/3) < 1 where the 
magnetization vanishes. 

The free energy in the ferromagnetic region can be computed at any necessary precision 
in a straightforward manner, it is sufficient in fact to expand Y^, a x ( a ) ^°s{ x { <T )) to the higher 
orders in A. In fact, at any order, the expansion is a compact analytic expression of A and 
n. 

Finally we remark that in the limit — > 0, 7 — > 00 and 7/3 = $ one recovers the fully 
connected mean field Ising model. In this limit, our expansion gives 6 C = 1, m c = A c , 
m c = 3(/3 — 1) and — (3F = |(/3 — l) 2 + log(2) which is the mean field Ising model solution 
in the same approximation. 



6 - DISCUSSION AND OUTLOOK 



In this paper we obtained the exact solution of the dilute Ising Model in the case of a 
vanishing temperature even in presence of a magnetic field. The solution gives a complete 
description of the geography of the associated random graph, not only the size of the con- 
tinent but also the number of isolated clusters of connected spins of all given size (islands) 
and the average number of bonds per spin both in the continent and in the islands. We 
were also able to find out the exact solution in the paramagnetic region and exactly deter- 
mine the phase transition line. On the contrary, in the positive temperature ferromagnetic 
region, our solution (13"T|) is the correct one only if the magnetization A c of the spins in the 
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continent is self-averaging in the thermodynamic limit. This assumption is not proven and 
the possibility that the magnetization has symmetry-broken characterization remains open, 
nevertheless, since at zero temperature there is not breaking and since the system is not 
frustrated, we are quite confident that our choice is correct. Furthermore, the ferromagnetic 
free energy for positive temperature is only given in terms of an expansion with respect to 
the order parameter A and, although the expansion can be made as accurate as necessary, 
the full and explicit free energy is still to be computed, if this is possible. 

Finally, we would like to mention that the present strategy can be straightforwardly 
applied to the dilute spin glass, but, unfortunately, it seems much more difficult to find 
analogous results. 
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